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ABSTRACT 

We review the available estimates of the masses of the compact object in Ultraluminous 
X-ray Sources (ULXs) and critically reconsider the stellar-mass versus intermediate- 
mass black hole interpretations. Black holes of several hundreds to thousands of M© are 
not required for the majority of ULXs, although they might be present in the handful 
of known hyper-luminous (~ 10^^ erg s~^) objects and/or some sources showing timing 
features in their power density spectra. At the same time, however, stellar mass BHs 
may be quite a reasonable explanation for ULXs below ~ lO"**^ erg s~^, but they need 
super-Eddington accretion and some suitable dependence of the beaming factor on 
the accretion rate in order to account for ULXs above this (isotropic) luminosity. We 
investigate in detail a 'third way' in which a proportion of ULXs contain « 30 — OOM© 
black holes formed in a low metallicity environment and accreting in a slightly critical 
regime and find that it can consistently account for the properties of bright ULXs. 
Surveys of ULX locations looking for a statistically meaningful relationship between 
ULX position, average luminosity and local metallicity will provide a definitive test of 
our proposal. 
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1 INTRODUCTION 

When, at the beginning of the '80s, point-like, off-nuclear X- 
ray sources were first detected in the field of nearby galaxies 
(see, e.g., |Fabbiano|[T989[ 12006] ), it was immediately recog- 
nised that the luminosity of a subset of these objects was un- 
usually large. If physically associated with their host galax- 
ies, these sources had an isotropic luminosity well in ex- 
cess of the Eddington limit for spherical accretion onto a 
IOMq compact object. These apparently Super-Eddington 
sources, later called UltraLuminous X-ray sources (ULXs), 



were first noticed in Einstein data (Long & van Speybroeck 
1983{ |Helfand|198"4l |M)biano 2006 ' . 



Nowadays well in excess of 150 candidate ULXs have 
been detected and catalogued by a variety of X-ray obser- 
vatories (see e.g. Roberts fc Warwick|2000 Colbert & Ptak 



2002 Swartz et al. 2004 Liu & Bregman 2005 1. A small frac- 



tion of ULXs are now known to be X-ray luminous interact- 



ing supernovae such as those described by Immler ( 2007 1 . A 



much larger fraction have subsequently been identified with 
background AGNs (~25%; , Swartz et al.|2004| see also|Fos-| 



chini et al. 2002b Masetti et al. 2003 Wong et al. 2008 1 



This background contamination is stronger in ellipticals, 
where it accounts for ~ 44% of all the ULXs, than in spi- 
rals (~14%; Swartz et al. 20041. However, the majority of 



these puzzling Super-Eddington sources constitute a very 
interesting class, that remains yet to be fully understood. 

The recent detection of a ~62 day modulation in the 
light curve of M 82 X-1 has been interpreted as the orbital 
period of the system ( Kaaret et al.|2006a|b Feng & Kaaret 



20071. A periodic modulation (12.5 hrs) has also recently 



been detected in another ULX (NGC 3379; Fabbiano et al, 
|2006j ). These results support the notion that ULXs are X- 
ray binary systems, where mass transferred from a donor 
star falls onto a compact object via an X-ray luminous ac- 
cretion disc. The high X-ray luminosities of ULXs suggest 
this compact object is most likely a black hole (BH). 

The X-ray spectral properties of ULXs show similar- 
ities with those of X-ray binaries (XRBs) in our Galaxy 
(e.g. [Foschini et al.||2002a[ ). In many cases the spectrum 
can be well reproduced by a multicolour disc (MCD) black- 
body, representing emission from an accretion disc, plus a 
power-law continuum (PL), with the latter nominally rep- 
resenting emission from a Compton-scattering corona. This 
is the same canonical model employed to describe the spec- 
tra of Galactic black hole X-ray binaries (cf. |McClintock| 
& Remillard 20061. Interestingly, the derived temperature 



of the MCD component in ULXs is often much lower than 
that observed in XRBs (e.g. [Miller et al.||2003[|2004b|[Feng| 
fc Kaaret]|2005 I . However, for some of the brightest ULXs, 
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a possible curvature above 2-3 keV has been reported and 
equally acceptable fits of their spectra may be obtained with 
(physically) diflferent models, that suggest the presence of 
hitherto unusual features such as an optically thick corona, 



a fast ionised outflow or a slim disc (e.g.ptobbart et al. 2006 
Gongalves Soria|2006l [Mizuno et al.|2007| ). 

Fast X-ray variability can also reveal much about 
accretion-powered sources. Although most ULXs show lit- 



tle variability on timescales of seconds to hours ( Swartz et 
|al.|2004| [Roberts et al.|2004[ ) , the recent detections of quasi 
periodic oscillations in the power density spectra of M 82 X- 
1 and NGC 5408 X-1 has shed new light on the timescales 



in the inner accretion disc of these systems ( 


Strohmayer & 


Mushotzky||2003| |Fiorito & Titarchuk||2004| 


Mucciarelli et 


al.||2006 


Strohmayer et al.|2007 


Casella et al.|2008 


I. 



Stellar optical counterparts have been discovered to be 
associated with a number of ULXs ([Roberts et al 



2001 



Goad et al.p002||Liu et al.|2 002, 2004; Kaare t et al.|2004a 



Soria et al. 



2005 



Zampieri eraL][2004l |Kaaret||200 5 Mucci arelli et al.|[2005 



PakuU et al. 



2006), although only some 
of them have been associated with stellar objects of known 
spectral type (e.g.pu et al.|2002l|2004| [Kaaret et al.|2004a[ 



Mucciarelli et al. 2005). In almost all cases, the counter 
parts appear to be hosted in young stellar environments 
(e.g. [Ramsey et al.|2006| [Pakull et al.|2006[ [Liu et al.|20"07l ) 

and have properties consistent with those of young, massive 
stars. However, some ULXs appear to be associated to older 
stellar populations and at least one possible later-type stellar 
counterpart is now known ( Feng fc Kaaret|2008 Roberts et 
al. 2008|. Some ULXs are also associated with very extended 
optical emission nebulae, that may provide important in- 
formation on the energetics and lifetime of these systems 
( Pakull fc Mirioni|2002[ [Roberts et al.|2003| ). These nebulae 
are also beginning to be detected as extended radio sources 
( Miller et al.|2005| [Lang et al.|200"7 l. 

The stellar environment of ULXs can also provide inter- 
esting constraints on the properties of ULX binary systems 
(ULXBs). Several ULXs are located in groups or clusters of 
OB stars. Isochrone fitting of the cluster colour-magnitude 
diagram has been attempted and provides cluster ages of 
tens of millions of years, although there is some disagree- 



ment among different authors (Ramsey et al. 2006 Pakull 
[et al.|2006[[Liu et al.|2007|[Gris6 et al.|2008[ ). These analyses 
translate into upper limits for the donor masses in ULXBs, 
assuming that they are coeval to their parent OB associa- 
tion. Typical values of these mass limits are in the range 
~ 10 — 2OM0 . Comparison of stellar evolutionary tracks of 
ULXs with the photometric properties of their optical coun- 
terparts on the colour-magnitude diagram may also be used 



to constrain the masses of their donor stars (e.g. Soria et al 
2005 Copperwheat et al.|2005 2007 1 . If accurate photome- 
try is available, this approach may also provide interesting 
clues to the BH mass, once binary evolution and X-ray irra- 



diation effects are taken into account (Patruno & Zampieri 



2008) 



ULXs play a fundamental role in the framework of X- 
ray source populations in nearby galaxies and can be de- 
tected and studied at larger distances than more 'normal' 
binary sources. An important tool to study the global prop- 
erties of these populations is their X-ray Luminosity Func- 
tion (XLF). [Grfmm et al.[ ( jSoOS] ) and [Gilfanov et al.[ ( |2004l ) 
found that the XLF of high mass X-ray binaries in the Milky 



Way, Magellanic Clouds and nearby starburst galaxies has 
a smooth, single power-law behaviour i n a broad luminos- 
ity range (10^^ — lO**^'^ erg s~^; see also Kaaret & Alonso- 
Herrero|2008 ). This suggests that ULXs with luminosity up 



to a few 10*" erg s~^ may simply represent the high lumi- 
nosity tail of the high mass X-ray binary population. 

These pieces of observational evidence, along with 
the long-term flux variability and the correlated luminos- 
ity/spectral variability (e.g. Kubota et al.[p()01| [La Parola[ 



et al. 2001 Zampieri et al. 



2004), strongly suggest that a 



large fraction of ULXs are accreting BH X-ray binaries with 
massive donors. The very high luminosity demands that the 
accretion rate be very high, even in case of efficient, disc 
accretion. A massive donor is then needed to fuel persistent 
ULXs, irrespectively of the BH mass (e.g. [Patruno et aT] 



20051 iRappaport et al. 120051 [Patruno & Zampieri|2008|), and 



the identification of blue, massive stars as the counterparts 
of some ULXs confirms this interpretation. However, tran- 
sient ULXs associated with older stellar populations, may be 
fueled through the rapid accretion of material accumulated 
in the accretion disc over a relatively long period of time, 
and not necessarily be associated with massive companions, 
as it is the case for, e.g, the Galactic BH candidate GRS 
1915-H05 ( [King|2002[ ). 

The critical issue is then understanding what is re- 
sponsible for the exceptionally high (isotropic) luminosity 
of these sources. Two main scenarios have been proposed. 
Firstly ULXs could be relatively normal stellar-mass BHs 
(< 2OM0) that are either anisotropically emitting X-ray bi- 



naries in a peculiar evolutionary stage (King et al. 2001 
King fc Pou"nds|[2003[ ), or are truly emitting above the Ed- 



dington limit via a massive, modified accretion disc structure 

two- 



(e.g. photon bubble dominated discs, Begelman|[2002 
phase super-Eddington, radiatively efficient discs, |Socrates| 
[fc Davis|[2006| slim discs, ,Ebisawa et al.^,2003[), or pe rhaps 



via some combination of the two ( Poutanen et al. 2007 King 



[2008 ). Secondly, the compact object could simply be big- 
ger, and the accretion would be in the usual sub-Eddington 
regime. In this case the compact object would be an inter- 
mediate mass black hole (IMBH) with a mass in excess of 
100 Mq (e.g. [Colbert fc Mushotzky[[1999[ ). Population syn- 
thesis calculations show that, in both scenarios, the mass 
transfer rates needed to supply the majority of the ULXs 
can be attained over a significant fraction of the life time 
of the systems and that the production efficiency of the two 
models are comparable (albeit with very large uncertainties 
on both), provided that stellar mass BHs can exceed the Ed- 



dington limit by a factor > 10 ( Podsiadlowski et al.[[2003 



Rappaport et al.|2005| [Madhusudhan et al. f2006, 20U8ir 



In this Paper we review the available estimates of the 
masses of the compact object in ULXs and present a critical 
re-evaluation of the current evidence regarding the stellar- 
mass versus intermediate-mass black hole interpretation. 



Here and in a companion investigation (Mapelli, Colpi fc 
[Zampieri[p009t we highlight an alternative formation sce- 
nario, already suggested before but never explored in depth, 
in which a proportion of ULXs contain 30 — 90Mq BHs 
formed in a low metallicity environment and accreting in a 
slightly critical regime. The plan of the paper refiects this 
approach. We start from a quite comprehensive review of 
the available estimates of the masses of the compact objects 
in ULXs (§ |2| and critically reconsider the 'traditional' in- 
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Table 1. Masses of the BH hosted in some ULXs estimated using X-ray spectroscopic methods. 





A, max 

^1040 ^iv..T fi — 1^ 
(iU erg s j 


Edd. limit'' 


MCD fif: 


Schwarzschild disc"* 


KERRBB fit^ 


SUm disc fit^ 








I AT \ 






M 81 X-1 


0.66 


130 






49+25 




M 81 X-9* 


1.1 


220 


330 - 470 


150 - 215 




M 101 X-2 


0.41 


82 






f;o + 127 

O'^-eo 




NGC 253 X-1 


0.29 


58 






73tfg 




NGC 253 X-3 


0.1 


20 






go + 113 




NGC 1313 X-1 
NGC 1313 X-2 
NGC 4559 X-7 
NGC 4559 X-10 


2.5 
1.5 
2.1 
1.2 


500 
300 
420 
240 


400l« 
2500tif„o 


200l« 

100+35 
-,200+1000 




16±1 
74±5 
311^^ 


NGC 5204 X-1 


0.5 


100 








23±3 



* Also known as Ho IX X-1 

" Cropper ct al. (2004); Kaarct et al. 72004a'); 'Mi ller et ah] | |2004a^ ; [Roberts et al.| | |2005| l; [Mizuno et al.| 
(|2007J ; ,Mucciarem et ah, ( ,2007J ; |Hui fc Krolik, 1 ,2008^ 
" ^BH computed from eqTTTT 

Mbh computed from eq. pf, wit h h = 9.5 and / = 1.7 ( [Miller et al.||2004at. The values of D and Kbb 
are taken from [Miller et al.|^004a[ ), [Cropper et al.[ ( [2b04[ ) and| Lorenzin fc Zampieri| (|2009|) (c ost = 1). 



A^BH computed from eq. ((2|, with h = 19.2 (Schwarzschild disc; Lorenziir^rZampieri"2009| and / 
The values of D and Kgg are from [Miller et al . ( 2004aj), [Cropper et al.[ ( [2004^ and Lorenzin 

(•2009) (cosi = 1). ' 

A/bh from the X-ray spectral fits of [Hui 



1.7. 
Zampieri[ 



Krolik|([2008t. 



^ A^BH from the X-ray spectral fits of [Vierd^anti^tal 
( [Vierdayanti et al.|2008^ . 



(20061, multiplied by a 1.3 correction factor 



terpretations of the nature of these sources (§|3|. We then 
discuss the low-metallicity scenario (§[4]) and some observa- 
tional tests to investigate it (§[5|. A conclusion section (§[6| 
follows. 



2 MASS ESTIMATES IN ULXS: METHODS 
AND RESULTS 

Thanks to the identification of the optical counterparts of 
a handful of ULXs, the measurement of the mass function 
of ULXBs is now well on the way to becoming possible, 
and will provide direct constraints on the masses of individ- 
ual sources. Unfortunately, given the observational difficul- 
ties associated with such measurements - most notably their 
faintness, with typical magnitudes in the range m„ ~ 22 — 26 
I Roberts et al.|2008 1 , and the contamination of the counter- 
part spectrum by nearby stars - only one very recent mea- 



surement of optical periodicity has been made (Liu, Breg- 



man fc McClintock||2009| ), and a mass function is yet to be 
constrained. However, this is where the observational effort 
is focussed at present and where the definitive answer to the 
question of whether intermediate or stellar mass BHs power 
ULXs will come from. Claims of a radial velocity shift of 
~ 300 km s~i in the He II A4686 line have been reported 
for the optical counterpart of NGC 1313 X-2 by [Pakull et] 
al. (2006), but this measurement may be uncertain (see e.g. 
MucciarelU et al.||2005| . 

Until these measurements are performed, we have to 
rely on indirect methods to estimate the BH mass. Assuming 
that the emission of ULXs originates from accretion and 
that it is stationary and isotropic, a lower limit for the BH 



mass Mbh is obtained for Eddington-limited accretion (L : 

^Edd): 



Ma 



10*0 erg s" 



200 



10*0 ergs" 



(1) 



where Lx typically refers to the [0.2-10] or [0.3-10] keV band 
and we consider a 'bolometric correction' of ~ 2 to account 
for the flux emitted outside this band. In fact, for the typ- 
ical spectral parameters of a bright ULX (colum density 
~ 3 X 10^^ cm~^ and power-law photon index ~ 1.7) the 
bolometric flux is ~ 2 times larger than the flux emitted in 
the 0.2-10 keV band (e.g. [Patruno fc Zamp icri 2008). This 
equation is also sensitive to the material accreted; the Ed- 
dington limit will rise by a factor ~ 2 for the accretion of 
helium (e.g. Grimm et al.|2003 |. Masses of some ULXs cal- 
culated from equation (|1| are reported in Table[l] Estimates 
based on this argument may become more reliable if there is 
some other evidence that the emission is (almost) isotropic, 
as for example when ULXs appear to be responsible for the 
photo-ionization of their surrounding optical nebulae (e.g. 
Kaaret et al.|2004a||Abolmasov et al.|[2()07| ). 

Spectroscopic estimates of the BH masses have also 
been attempted assuming that the soft component observed 
in ULX X-ray spectra can be ascribed to emission from 
an accretion disc (e.g. Miller et al. 2003[ [2004a). Assum- 
ing that this spectral component can be modelled with the 



so-called multicolour disc blackbody model (MCD; Mitsuda 
et al.||1984 l, it is possible to express the BH mass Mbh as 



(e.g. Lorenzin fc Zampieri|2009 1 



Mb 



,67.5 



Mp 



D 



IMpc 



fKBB\ 

\ cosi / 



1/2 



(2) 



where D is the distance of the source, / is a colour correc- 
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Table 2. Masses of the BH hosted in some ULXs estimated using timing methods. 





A ,Tnax 

(lO^O erg s-i) 


Edd. limit'' 


QPO" 
(Mq) 


Break/Lack of variability'* 
(Mq) 


Ho II X-1 


1.7 


340 




< 100 


M 82 X-1 


17 


3400 


95-1300 


25-520 


NGC 5408 X-1 


0.85 


170 


115-1300 


~ 100 


NGC 4559 X-7 


2.1 


420 




38-1300 



" From |Cropper et al.| f2004t ; [Kaaret et al.| | (2004a| ; iMucciareUi et al.| l |2006| ; [Strohmayer et al.| | |2007[ |. 
^ Mbh computed from eq. ||l||. 



Casella et al.| I 


2008|l. 


Cropper et al. 


(|2004|l 



tion factor that accounts for transfer effects (e.g. Shimura 
fc Takahara||T995| [Zampieri eraL]|200T| [Davis et al.||2005[ 
Hui et al.|2005 1 , 6 is the inner radius in units of the gravita- 
tional radius, and Kbb is the MCD normalization inferred 



from the spectral fit. In their work Miller et al. (20031 and 
[Miller et al] ( [2004a[ ) adopted b = 9.5 and / = 1.7 and de- 
rived estimates of Mbh well in excess of IOOMq for M 81 
X-9, NGC 1313 X-1 and X-2. A similarly large value of Mbh 
is obtained also for NGC 4559 X-7 using the MCD spectral 
parameters of Cropper et al. (20041 (see Table [l]). However, 
recently Lorenzin & Zampieri ( 2009 1 computed appropriate 



correction factors for b in the case of a relativistic, standard 
accretion disc and showed that, unless the BH is maximally 
rotating, the BH masses inferred for the same sources can 
be significantly lower than the values estimated by [Miller[ 
et al.[ ( |2003[ ) and [MiUer et ah] (|2004a[) (see again TakejlJ). 
Smaller BH masses have also been obtained through direct 
spectral fits of relativistic accretion disc models for a sample 



of disc-dominated ULXs by Hui & Krolik ( 2008 1 



The interpretation of the soft component in terms of 
emission from a standard accretion disc suffers from a high 
degree of equivocality. Spectral fits with a disc component 
and a comptonizing corona to the best available X-ray spec- 
tral data indicate that ULXs display distinct spectral cur- 



vature above 2 keV (Stobbart et al. 2006 Gladstone et al 



2009). This is because the corona is optically thick and cool. 



and hence hides the inner part of the accretion disc, in what 
is likely an extreme form of the so-called very high state 
of Galactic BH candidates ( Done fc Kubota|2006 Stobbart 



[et al.[[2b06[ ). In this case, the soft component is produced 
by the visible outer regions of the accretion disc. Only up- 
per bounds to the inner disc radius can be obtained and, 
similarly, the spectroscopic estimates for Afsn reported in 
Table [l] should be considered as upper limits. However, as 
the thick corona is probably only launched at extreme ac- 
cretion rates, this implies that this "ultraluminous" state is 
in the Super-Eddington regime, and hence the BH masses 
are r elatively smaU (< IOOM q; |Roberts[[2007) . Further to 
this, [Gladstone et al.| ( [2009[ ) consider the energy required 
to launch the thick coronae, and from that calculate the 
intrinsic (corona-less) disc temperatures, mainly recovering 
temperatures in the correct regime for the discs around 
stellar-mass BHs (0.7 - 1 keV). A subset of ULXs retain 
apparently cool discs even after this correction; however, 
[Gladstone et al.[ ( |2009| ) argue this is because these are the 
highest accretion rate stellar-mass BHs, in which a strong 



wind is launched from the central regions, creating a cool 
photosphere. In a similar vein, for other models character- 
ising super-Eddington accretion, such as the slim disc (e.g. 

Ebisawa et al.[[2003 1 or photon bubble 



Watarai et al. 2000 



models \Fmke fc Bottcher|2007[ ), the entire 0.2-10 keV spec- 
trum is produced in the accretion disc, although its physical 
state is completely different with respect to that of a stan- 
dard disc. However, indirect estimates of Mbh can again be 
obtained from X-ray spectral fits and give values typically 
in the range ^ 15 ~ 75Mq ( [Vierdayanti et al.|[2006l [2008 
see again Table [l]). 

The mass estimates inferred from X-ray spectral fits of 
ULXs depend critically on the interpretation of their spectra 
which, as mentioned above, is not unique. The situation will 
improve in the future as our understanding of the spectral 
evolution of ULXs will increase and it will be possible to 
select spectral models on the basis of their consistency with 
the observed correlation patterns (see e.g. [Feng fc Kaaret[ 
2009; Kaj ava fc Poutanen|2009] for preliminary work on ULX 
spectral variability based on XMM-Newton data). Nonethe- 
less, the available data points towards BH masses definitely 
smaller than those estimated from the early MCD spectral 
fits. 

In the last few years. X-ray timing has provided a new 
opportunity to estimate BH masses in ULXs thanks in par- 
ticular to the detection of broad band noise and quasi pe- 
riodic oscillations (QPOs) in the power density spectrum 
(PDS) of some ULXs, as M 82 X-1 (vqpq = 54 - 166 
mHz; [Strohmayer fc Mushotzkyj [2003} [Mucciarelli et al 
2006]) and NGC 5408 X-1 {uqpo ^ 20 mHz; [Strohmayer 
et al.[[2007 l. Extrapolating timing and spectral-timing cor- 
relations known to exist for similar timing features in BH 
binaries and assuming that the frequency of the QPO scales 
inversely to Mbh, various estimates have been obtained, 
which are consistent with a rather broad interval of values: 
10- IOOOMq for M 82 X-1 ([Fiorito fc Titarchuk|2004[ puc^ 



ciareUi et al.| 2006|[Feng fc kaaret|2007[ ), 600-5000 for NGC 
5408 X-1 (Str ohrnayer et al.[[2007[ ). Recently, a new timing 
approach has been presented to assess BH masses in ULXs, 
which is based the so called "variability plane", populated 
by both Galactic black-hole candidates and active galactic 
nuclei. Assuming that the accretion flow in ULXs behaves in 
a similar way (which remains an open question) and taking 
into account the uncertainty on the efficiency of the accre- 
tion disc, [Casella et al.| ( [2008[ ) find that Mbh is in the inter- 
val ~ 95 - 1300Mq for M 82 X-1 and ~ 115 - 13OOM0 for 
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NGC 5408 X-1 (see Table [2]|. In combination with QPOs, 
the scaling of the break frequencies of the broad band noise, 
by comparison with the corresponding timing features of 
Galactic BH candidates, has also been proposed for esti- 
mating Mbh- This method has been applied to NGC 5408 
X-1 and gives a similar range of masses to QPOs (see Ta- 
ble [2|. A tentative identification of a break at ~ 28 mHz in 
the PDS of NGC 4559 X-7 has also been repo rted ( [Cropper 



et al.||2004| , although [Barnard et aL] ( |2007| ) called it into 



question. The inferred BH mass is reported in Table [2| 

The problem with using timing properties of ULXs to 
directly infer masses is that it remains to date unclear how 
exactly these quantities are related. All the estimates are 
based upon tentative identifications of the timing features 
and the use of scaling laws that are known to hold only 
for a limited number of objects, and are therefore highly 
uncertain. 

A further possibility is to use the non-detection of vari- 
ability power in the PDS to limit the size of the BH, as- 
suming that all power is at higher frequencies, as is seen 
in various XRB states. Adopting this approach [Goad et[ 



all (|2006| found Mbh < IOOMq for H olmberg H X-1 (see 
again iablej^l. In fact, Heil et al. (20091 have demonstrated 
that suppressed temporal variability (compared to the PDS 
of classic XRBs and AGNs) appears a common feature of 
ULXs. There are several possible explanations of this - the 
variability is limited to higher frequencies, the data in the 
XMM-Newton band pass are disc-dominated (and so any 
variability is heavily diluted) or the ULXs are in a new, 
super-Eddington accretion state in which the X-ray emis- 
sion is very stable (note this is predicted in the hydrody- 
namic simulations of highly super-Eddington accretion by 
Ohsuga[[2007[ ). A gain, the common thread running through 
all these models is that the BH is relatively small. 



3 INTERMEDIATE OR STELLAR-MASS BHS? 

We have reviewed the mass estimates drawn from obser- 
vations of individual ULXs; we now ask how these results 
fit into our more general understanding of the possible na- 
ture of their underlying engines. ULX models differ mainly 
in the assumptions on the physical state of the disc and 
its mode of accretion. If the accretion disc is in a standard 
regime, then emission is isotropic and the most straight- 
forward interpretation of the exceptionally high luminosity 
of ULXs is that they contain BHs with large masses. Both 
the very high luminosity and low characteristic temperature 
(and high normalization, see previous Section) of the soft 
spectral component have been taken as evidence for this in- 
terpretation. 

But how big is the BH mass? Early estimates based 
on equations ([T]) and ([2| gave masses largely in excess of 
IOOM0 (up to several thousands; see Table [T]). The obvious 
question is then how a BH this massive may form. It has 
been proposed that remnants from the collapse of Popula- 



tion III stars formed in cosmological epochs ( Madau & Rees 
[2001[ ) may trigger ULX activity, if they can capture a donor 
star to accrete from. A second formation route for IMBHs 
may be in globular clusters, through repeated mergers of 
stellar mass BHs ( [Miller fc Hamilton[[2002| ), or m young, 
dense stellar super clusters, from the dynamical collapse of 



supermassive stars in their centres (e.g. Portegies Zwart et 
al.[[2004 1. ULX activity would be sustained by binary com- 



panions captured in the cluster (Blecha et al. 20061. The 



latter may be a possible explanation for some ULXs (e.g. 
M82 X-1) but, in galaxies with starburst activity, the ma- 
jority do not appear inside such supermassive clusters (e.g. 



Zezas et al. 2002; Kaaret et al. 2004b I. Another difliculty 



with the IMBH interpretation is the apparent break in the 
power-law slope of the XLF of the high mass X-ray binary 
population in external galaxies at a luminosity ~ 2x 10*" erg 
s~^. If there is a population of large IMBHs contributing to 
the XLF, the break suggests that they turn off at this lumi- 
nosity; yet this is at a rather low fraction of the Eddington 
limit for putative large IMBHs (~ 10% for a IOOOMq BH). 
This would be rather unusual behaviour, as no other accret- 
ing class switches off at only a fraction of their Eddington 
limit ( Roberts|2007 |. Finally, the co-location of ULXs with 
regions of star formation, such as those in the Antennae and 
Cartwheel galaxies, implies that they must be (relatively) 
short-lived, which requires successive generations of ULXs 
to be formed over the duration of the star formation event. 
Thus, if all ULXs in these regions were IMBHs, an unfea- 
sibly large fraction of star forming mass would end up in 
IMBHs ( [King|2004| [MapeUi et al.|2008[ ). In principle, these 
arguments rule out all but a small minority of ULXs from 
being IMBHs bigger than ~ IOOMq. 

So can we explain ULXs as stellar-mass BHs? If the 
accretion flow in ULXs is in a different regime, the situ- 
ation may be different and either the isotropy and/or the 
Eddington limit may be circumvented. This occurs if the 
accretion rate is at or above Eddington, so that radial advec- 



tion of thermal and radiative energy (slim disc; Abramow- 



icz et al. 1988 Ebisawa et al. 20031 and/or radiation- 



driven instabilities (photon bubble model; [Begelman|[2002| 



2006 1 set in. A slim disc can sustain larger accretion rates 



than a standard disc and modest super-Eddington lumi- 
nosities (< lOI/Edd)- As the emission is isotropic, luminos- 
ity scaling arguments similar to those discussed above give 
A/bh > 2O(Lx/lO*°ergs"^)M0. Therefore, the presence of 
slim discs would imply that only the brightest known ULXs 
(with Lx ^ 3 X lO*'' erg s^^) need BHs significantly more 
massive than the stellar-mass BHs in our Galaxy (see also 
Table[TJ . Accretion discs dominated by photon bubble trans- 
port may also reach super-Eddington luminosities, while re- 
maining geometrically thin. For a 2QMq BH, the maximum 
luminosity is ~ 30LBdd ( [Begelman|2006| ) . This may in prin- 
ciple account for the emission of all but the very brightest 
(> 5 X 10*° erg s~^) ULXs, but photon bubble-dominated 
discs are subject to the same thermal and viscous instabili- 
ties that characterize the inner region of radiation pressure 
dominated discs and hence may be significantly unsteady 
(with more than a factor of 10 variation in the emitted flux) 
on short time scales (e.g. Zampieri et al. 20011. However, 
bright (> 10*° erg s~^) ULXs typically show more limited 
fluctuations in the observed X-ray luminosity. An alterna- 
tive scenario for generating steady, super-Eddington lumi- 
nosities from stellar-mass BHs, is the radiatively efficient, 
two-phase super-Eddington accretion disc model by |Socrates| 



& Davis ( 2006 1. In this model the gravitational potential en- 



ergy is not trapped in the disc, but effectively removed from 
it through magnetic buoyancy and dissipated in a corona. 
Fields anchored in the disc and/or Compton drag in the 
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low density corona prevent the launching of a wind, keeping 
the radiative efficiency high and assuring super-Eddington 
luminosities. However, we note that there are many assump- 
tions and theoretical uncertainties that need to be clarified 
in this model, and a quantitative estimate of the maximum 
attainable luminosity is not yet available. 

If accretion becomes largely super-Eddington, other 
processes may be present that complicate the picture. A 
thick disc may form and the emission becomes beamed 



(King et al. 2001 King 20021. At the same time, outflows 



and powerful ejection of matter along the axis perpendic- 
ular to the accretion disc may be produced (as in SS433; 
e.g. Poutanen et al. 20071. In these assumptions, emission 
is no longer isotropic. If Li^o is the apparent isotropic lu- 
minosity, the BH mass inferred from equation |l| must 
then be corrected for the so called beaming factor bf, 
that represents the fractional opening angle of the beam: 
Mbh =^ 2O(6//O.l)(Lo.2-io,iso/lO*°ergs-i)M0. Hence, sim- 
ple beaming can not account for bright (> 2 x 10*° erg 
s~^) ULXs, unless one is willing to consider a beaming 
factor < 0.05 and therefore a very narrow opening an- 
gle (< 18°), that appears to be more consistent with a 
jet rather than a geometrical funnel in a thick disc. How- 
ever, in addition to having beamed emission, thick discs 
may also radiate at super-Eddington lur ninosities, reach- 



ing at most (1 -I- ^n(M/MEdd))iEdd (e.g. Poutanen et al. 
20071 |King]|2008[ ). A hypothetical ULX in this state would 
bear some similarity to a Galactic BH candidate in the very 
high state, albeit with the ULX being in a much more ex- 
treme version of this state, with powerful winds carrying 
away the excess matter and energy, potentially thickening 
the corona and even producing a cool photosphere (consis- 
tent with the X-ray spectral modelling of [Gladstone et al.| 
2009). A combination of a beaming factor bf ~ 0.3 — 0.5 
and super-Eddington emission (~ WLEdd) may in principle 
explain ULXs with luminosities up to « 10*° erg s~^, as- 
suming accretion onto a stellar-mass BH. Furthermore, if a 
dependence of the beaming f actor on the accretion rate is 
assumed (6/ oc (M/A^Edd)~^; King 2009), it might be pos- 
sible to account also for the high luminosity tail of the ULX 
population. 

In order to assess the viability of the different scenar- 
ios for the origin of ULXs, it is necessary to understand the 
possible evolutionary history of the various types of candi- 
date binary systems and compare them with the available 
observations. Galculations of the evolutionary tracks of ULX 
binaries and model population studies of systems contain- 
ing stellar-mass BHs and IMBHs have been carried out by 
several authors (e.g. [Podsiadlowski et al.||2003| [Patruno et| 



aL]|2005[ [Rappaport et al.||2005| [Madhusudhan et al.|2006 ) 
Although calculations depend sensitively on uncertain pa- 
rameters of the common envelope phase, it turns out that 
stellar-mass BHs accreting at super-Eddington rates may 
be able to account for most of the observed ULXs (except 
for the brightest), if they violate the Eddington limit by a 
factor ~ 10 — 30 ( Podsiadlowski et al.||2003 Rappaport et 



al. 2005). Similarly, IMBH systems might produce bright 
[Lx ^ 10*° erg s~^), persistent ULXs and have an accept- 
able production efficiency if the donor star is > IOMq and 
the initial orbital separation is small (< 6 — 40 times the ini- 
tial main sequence radius of the donor; [Patruno et al.|2005l 
Madhusudhan et al.|2006l ). 



Along the same lines, theoretical calculations of the 
color-magnitude (CM) diagrams for systems containing 
stellar-mass BHs and IMBHs are used to constrain ULX 
models and their evolutionary history by comparison with 



observations of their optical counterparts ( Madhusudhan et 
al.|2008 Patruno Zampieri|2008 ). As already mentioned, 
in order to supply the mass transfer rates needed to fuel 
ULXs, rather massive donor stars are required. The evolu- 
tionary tracks of such systems are strongly affected by the 
binary interaction and the emission from the accretion disk, 
including X-ray irradiation. Numerical computations show 
that the regions of the CM diagram with the highest prob- 
ability of finding ULX optical counterparts have B — V be- 
tween ~ —0.1 and —0.3 and correspond to the early phases 
of the evolution of massive donors, while they are on the 
main sequence or the subgiant branch ( Madhusudhan et al. 
^2008) . This result is consistent with the properties of the ob- 
served counterparts. Similarly, the most favourable orbital 
periods are between 1 and 10 days, corresponding again to 
the main sequence or subgiant phases. This is true for both 
stellar-mass BH and IMBH systems and depends on the fact 
that the donors spend most of their life time in these phases. 
In these conditions normal nuclear-driven mass transfer is ef- 
fective and provides sufficiently high transfer rates to sustain 
the ULX emission, although in rare circumstances the evo- 
lution may be driven by the thermal timescale mass transfer 
during the giant phase. 



4 A DIFFERENT INTERPRETATION 

A critical revaluation of the available observational evidence 
presented in Sections [2] and |3] indicates that BHs of several 
hundreds to thousands Mq are not required for the majority 
of ULXs. However, it is not possible to rule out that they 
are present in the handful of known hyper-luminous (~ 10*^ 
erg s~^) objects and/or in the sources showing large ampli- 
tude broad band noise or QPOs in their PDS (such as M 
82 X-1 and NGC 5408 X-1; see Table[2|. At the same time, 
models with stellar mass BHs may work for a large fraction 
of the ULX population, if the accretion flow has some de- 
greee of beaming and is super-Eddington. Bright (> 10*° erg 
s~^) ULXs may be accounted for if some form of modified 
beaming that accounts for a suitable dependence of fe/ on 
the parameters of the accretion fiow is allowed ( |King|20"09 ) . 
Although none of these scenarios can be ruled out, the fact 
that the observational limits discussed in Section |2] are con- 
verging towards masses < IOOM0, but bigger than stellar 
mass BHs, led us to consider an alternative interpretation. 
In our scenario bright ULXs may contain BHs with masses 
above 30-40 Mq and up to ~ 80 — 9OM0 , formed from ordi- 
nary stellar evolution of massive (30 — I2OM0) stars in a low 
metallicity natal environment. While this idea has already 



been suggested before (e.g. PakuU & Mirioni 2002 Crop- 
per et al. 2004 Zampieri et al. 2004 1 , it has not yet been 



explored quantitatively in detail. 

Stars with mass > 8Mq produce compact remnants 
from the gravitational collapse of the iron core. For stars up 
to ~ 25 — 3OM0 the collapse is halted when the core reaches 
nuclear densities: the star explodes and a neutron star forms. 
For larger main sequence masses, the early acccretion of 
the inner mantle onto the core before shock passage and 
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Main sequence mass (Mq) 



Figure 1. Final mass as a function of initial mass for stars of different metallicities as computed by Maeder (1990, 1992; black), Maeder 
& Myenet (2001; blue), Portinari et al. (1998, red), and Eldridge & Tout (2004, cyan). 



the fallback of material afterwards cause the newly formed 
proto-neutron star to collapse to a BH after the supernova 
explosion ( |Zampieri||2002| and references therein) . At solar 
metallicity, these fallback BHs reach at most ~ IOMq as, 
for very massive stars with mass > 40Mq, the stellar enve- 
lope is in large part effectively removed through line-driven 
winds, while the remaining part is expelled by the supernova 
explosion. 

For sub-solar metallicities, however, this mechanism 
becomes progressively less efficient and stars with masses 
above ~ 30 — AOMq may retain rather massive envelopes at 
the time of explosion. The supernova shock wave then loses 
more and more energy in trying to unbind the envelope until 
it stalls and most of the star collapses to form a BH with 



a mass comparable to that of the pre-supernova star ( Fryer 



1999). These may be the BHs hosted in some ULXs. Their 
mass would not be significantly larger than ~ 80 — 9OM0 
as above ~ 100 — I2OM0 a star undergoes pulsational pair- 
instability in its core and most of the envelope mass is ex- 
pelled. We note that the possibility of forming BHs in this 
mass range through a different channel (binary mergers of 
massive components) was also proposed few years ago (e.g. 
Belczynski, Sad owski fc Rasio|2004[ ). 

Computation of the evolution of massive stars up to ad- 
vanced evolutionary stages for different metallicities and/or 
including mass loss have been performed by several authors 



(e.g. 


Hellings & Vanbeverenll981 


■Macder'1990 


1992 Porti- 


nari et al.|1998||Heger et al.|2003 


,Chicffi and Limongi|2004| 


Eldridge and Tout 2004 


Hirschi|2007 1. These works adopted 



known empirical parameterizations of the mass loss rate for 
stars over the whole Hertzsprung-Russell diagram (e.g. [de| 
Jager et al.||1988 Nieuwenhuijzen fc de Jager||1990 1. A cer- 



tain degree of uncertainty is introduced if the star enters 
some peculiar evolutionary stages, such as the Wolf-Rayet 
(WR) stage (e.g. |Langer|[T989l [Wellstein fc Langer|p^ . 
In fact, mass loss rates in WR and O stars may be signif- 
icantly reduced if the wind is clumpy as a consequence of, 
e.g., supersonic turbulence ( [Moffat fc Carme llc 1994; Fuller- 1 
ton et ar]|2006 1. In particular, during the WR phase a de- 



crease by a factor of > 3 with respect to homogeneous wind 
models is attainable. This turns out to be compatible with 
some observational estimates and would clearly lead to more 
massive pre-supernova stars and hence more massive BHs. 
Additional uncertainty is caused by the dependence of mass 
loss on metallicity. A scaling law oc Z°-^ is often adopted 
for hot stars (see e.g. Kudritzki et al.|19 89, Nugis fc Lamers 
[21300 j ). For example, at the end of main sequence, the mass 
of a star with an initial mass of IOOA^q is ~ 25Af0 at solar 
metallicity and ~ 75 Mq for Z ~ O.IZq. However, the mass 
lost during the He burning phase may be much more signif- 
icant than that occurring during main sequence. According 
to the adopted mass loss history, the final mass of a star 
(after all nuclear burning stages are exhausted) may differ 
up to a factor of ~ 2 (or even more for clumpy winds) for 
a given metallicity. Considering again a star with an initial 
mass of IOOM0, Figure [1] shows that its final mass may be 
in the interval ~ 3 - 6Mq for Z ^ Zq and ~ 30 - 7OM0 for 
Z^O.IZq. 

As already mentioned above, the final evolutionary 
stages of a star and, in particular, the outcome of the final 
collapse depend critically on how massive is the envelope 
that it retains at the time of explosion. Therefore, owing 
to their larger final masses, the fate of stars with sub-solar 
metallicity is likely to be be quite different from that of 
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higher metalhcity stars. Ahhough different authors obtain 
different results for the mass of the compact remnant, it is 
not unreasonable to think that, if an envelope more massive 
than ~ 30 — AOMq is retained at the time of explosion, a low 
metallicity {Z ~ O.IZq) star may collapse directly to form a 
BH of comparable mass. Looking again at Figurejl] it is pos- 
sible to see that, already at metallicities < 0.2Zq, stars with 
initial mass > 6OM0 appear to possess final envelope masses 
that overcome this threshold for direct BH formation. Sig- 
nificant stellar rotation (hundreds of km s^^) may change 
this picture somewhat, as it enhances the mixing of heavy 
elements throughout the star increasing the metal content 



of the envelope and consequently mass loss (e.g. Maeder & 
Meynet|200l||Meynet fc Maeder|2005| ). Also, ejection of part 
of the envelope during the final collapse can not be ruled out. 
However, if the core is not rapidly rotating, there is no good 
reason why most of the star should not collapse into a BH. 
Therefore, the formation of a > 30 - 4QMq BH through- 
out the evolution of a low metallicity, slowly rotating star 
of 40 — 120Mq appears a viable possibility. The parame- 
ter space in the metallicity-main sequence mass plane where 
this formation channel may actually work corresponds to the 
black-shaded area between ~ 40 and ~ lOOM© in Figure 1 



2004 



3gei 

Er 



of Heger et al. ( 2003 1 (see also Figure 5 in Eldridge and Tout 



At variance with intermediate mass BHs, the formation 
of these very massive stellar remnant BHs does not require 
an exotic, new mechanism but is referable to ordinary stellar 
evolution. Also, the unbroken power-law slope of the X-ray 
binary population up to ~ 2 x 10*" erg s~^ is consistent 
with the fact that variations in metallicity may produce a 
continuum distribution of BH masses from the stellar-mass 
BHs of 10 - 2OM0 up to the suggested BHs oi^40- QOMq . 
Given the size of these BHs, no difficulty with the fraction 
of star-forming mass in large starbursts ending up in BHs 
would arise, and so the objection of |King| (j2004) for ULXs as 
~ IOOOMq IMBHs is circumvented. At the same time, only 
modest beaming (bf ~ 0.5) or slight violations of the Ed- 
dington limit (a factor of a few) would be needed to account 
for the luminosity of bright (> 10*" erg s~^) ULXs, at vari- 
ance with the extreme accretion scenarios required by stellar 
mass BH models. Also the essentially isotropic irradiation of 
X-ray photoionised nebulae would find an explanation. 



5 OBSERVATIONAL TESTS 

Model population studies for our scenario are needed to 
determine the production efficiency of binary systems con- 
taining very massive BHs and to understand if they are in 
agreement, in a statistical sense, with the available X-ray 
and optical data of ULXs. However, there may be already 
some indications strengthening our present suggestion. In 
a parallel, preliminary investigation, we show that massive 
BHs formed in low metallicity environments might well ex- 
plain most of the ULXs observed in the Cartwheel galaxy 



(Mapelli, Colpi & Zampieri 20091. Also, the optical lumi- 



nosities of massive BH systems would be, on average, larger 

^ An estimate of the remnant mass from massive stars and its 
dependence on metallicity and wind mass loss rates has been very 



than that of stellar-mass BHs, as the former allow for more 
massive donors (> 25Mq) and have more extended accre- 
tion discs that dominate the optical emission. This would 
make them more consistent than stellar-mass BHs with the 
observed distribution of the luminosity of the ULX opti- 
cal counterparts (e.g. Madhusudhan et al.|2008 Patruno & 
Zampieri|20 08| 



A crucial aspect of the interpretation of ULXs in terms 
of BHs from the direct collapse of \ow-Z, massive stars is 
the metallicity of the environment in which ULXBs form. 
The available estimates appear to favour a low metallicity 
scenario, although there are some discrepancies. Optical ob- 
servations appear to provide evidence of sub-solar metallic- 
ity in the environment of some ULXs. The emission neb- 
ula surrounding Ho H X-1 has a spectrum resembling that 
of a high-excitation HII region typical for low metallicity 
{Z ~ O.IZq) star-forming regions ( |Mirioni|[2002l |Pakull k 
|Mirioni|[20^ . The stellar environment of NGC 4559 X-7 
shows a blue-to-red supergiant ratio (3 ± 1) and colours of 
the red supergiant population consistent with a low metal 
abundance environment with Z = 0.1 — OAZq, similar to 



that in the SMC and other nearby dwarf galaxies ( Soria et 



alj[2005j). Also the stellar field around NGC 1313 X-2 is 
characterised by a low metallicity {[Fe/H] = —1.9 ± 0.3), 
as inferred from the intrinsic colour of the red giant branch 
(Gri se et al.||200^ l; studies of the metal abundance of HII 
regions in NGC 1313 also give low values of Z (~ 0.008; 
Walsh fc Roy|p97l [Hadfield fc Crowther|2007 1. 



Winter et al 



( |2007| l analysed high signal-to-noise XMM 



recently derived also by Belczynski et al. 1 2009 I 



spectra of a sample of 14 ULXs, trying to determine the Oxy- 
gen abundance from the detection of K-shell photoioniza- 
tion edges. They apparently find values that match the solar 
abundance. However, the comparison of the X-ray estimates 
with a compilation of [O/H] ratios determined through spec- 
trophotometric studies of HII regions ( [Pilyugin^ 2004) and 
with the luminosity-metallicity relation derived from the 
Sloan Digital Sky Survey ( Tremonti|2004 1 shows significant 
systematic differences, especially at low galaxy luminosities 
and sub-solar metallicities. We note that Ho II X-1 provides 
a clear example of this dichotomy for a ULX; its XMM- 
Newton RGS spectrum suggests a significantly higher metal- 
licity (~ 0.6 times solar; Goad et al.|2006" l than the optical 
data. At the same time, X-ray spectral fits in at least the 
case of NGC 4559 X-7 seem to provide evidence for a more 
subsolar metallicity (Z ~ 0.3^0 ; Cropper et al.|2004 |. 

It will be possible also to test our proposal against the 
stellar-mass BH interpretation as it will lead to a differ- 
ent spatial distribution of bright ULXs. In fact, ULXs from 
stellar-mass BHs should essentially appear anywhere in re- 
gions of star formation or in young stellar environments, 
regardless of metallicity. Indeed there may even be a bias 
towards these objects appearing in low metallicity regions, 
as effects such as the low mass loss rate in stellar winds will 
tend to keep binaries close, allowing more high mass trans- 
fer systems. However, if we make the reasonable assumption 
that the modes of accretion and Eddington rates will be 
similar in both standard ~ 10 — 2OM0 black holes, and the 
~ 30 — 9OM0 black holes examined in this paper, then the 
latter black holes may simply be on average brighter. Hence, 
in our proposed scenario ULXs should show some evidence 
of correlation (in terms of position and average luminosity) 
with low metallicity environments. So, one of the definitive 
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tests of our proposal would be to survey ULX locations, and 
determine whether a relationship between ULX luminosity 
and local metallicity was evident in a large enough sample 
to provide statistically meaningful results. 

We note that evidence supporting our argument is al- 
ready available. Firstly, Swartz et al. (20081 have recently 



surveyed galaxies within the Local Volume and determined 
that the specific ULX frequency decreases with host galaxy 
mass above ~ IO^ '^Mq. This means that smaller, lower 
metallicity systems have more ULXs per unit mass than 
larger galaxies, consistent with the idea that BHs can at least 
form and/or feed more efficiently in low metallicity environ- 
ments. Secondly, we note that there is interesting evidence 
from our own relative backyard that the brighter ULXs are 
more abundant in the late-type spiral galaxies we might ex- 
pect to be low-metallicity systems. The ROSAT High Res- 



olution Imager ULX survey of Liu & Bregman ( 2005 1 lists 
15 ULXs within 5 Mpc|^ four of which have observed X-ray 
luminosities in excess of 5 x 10^^ erg s"'^. All four of the very 
luminous ULXs reside in galaxies of Hubble type Sd or later. 
This compares to only 3 of the 11 lower-luminosity ULXs be- 
ing hosted by similarly late systems - the remainder are in 
galaxies of type between Sab - Scd. Although the sample of 
ULXs is small, this clearly supports the case that brighter 
ULXs preferentially occur in the smaller, lower- metallicity 
systems where we might expect to find the very massive 
stellar remnant BHs. 



It is worth noting that, recently, Prestwich et al. (20071 



and Silverman & Filippenko ( 2008 1 succeeded in perform- 



ing dynamical mass measurements using Gemini and Keck 
spectra of the Wolf-Rayet optical counterpart of IC 10 X-1, 
a variable X-ray source in the the Local Group metal poor 
starbust galaxy IC 10. They find a BH mass in the range 
23 — 33Mq, which represents the most massive BH known 
to exist in a binary system and definitely corroborates our 
interpretation. 



6 CONCLUSIONS 

In the last few years. X-ray and optical observations have 
significantly boosted our understanding of ULXs. We are 
now confident that the majority of these sources are X-ray 
binaries in external galaxies and we suspect that many may 
have massive binary companions. Yet, the most fundamental 
questions on ULXs still remain to be definitively answered: 
do they contain stellar or intermediate mass BHs? How do 
they form? 

A critical revaluation of the available evidence presented 
indicates that BHs of several hundreds to thousands Mq 
are not required for the majority of ULXs, although they 
might be present in the handful of known hyper-luminous 
(~ 10*^ erg s~^) objects and/or some sources showing tim- 
ing features in their power density spectra. At the same 
time, however, stellar mass BHs may be quite a reasonable 
explanation for ULXs below ~ 10"*" erg s~^, but they need 
super-Eddington accretion and some suitable dependence of 
the beaming factor on the accretion rate in order to account 
for ULXs above this (isotropic) luminosity. 



^ Their survey covers 27 galaxies within 5 Mpc, 17 of which are 
of Hubble type Sd or later. 



We investigated in detail an alternative scenario in 
which bright ULXs contain BHs with masses above ~ 30 — 
AQMq and up to ~ 80— 9OM0, produced by stars with initial, 
main sequence mass above ~ 40 — SOMq. At sub-solar metal- 
licity, the explosion energy of these stars is not sufficient to 
unbind the envelope and most of the star collapses to form 
a BH with a mass comparable to that of the pre-supernova 
star. These may be the BHs hosted in bright ULXs. Above 
~ 100— 120Mq pulsational instability becomes effective and 
most of the envelope is expelled from the star. 

The formation of these very massive stellar remnant 
BHs does not require an exotic, new mechanism but is refer- 
able to ordinary stellar evolution. For luminosities ~ 10*'^ 
erg s~^, this would imply only modest violations of the 
Eddington limit, attainable through very modest beaming 
ipi ~ 0.5) and/or slightly super-critical accretion. 

Measurements of the metallicity of the environment of 
some ULXs appear to favour a low metallicity scenario, al- 
though there are some discrepancies. Surveys of ULX loca- 
tions looking for a statistically meaningful relationship be- 
tween between ULX number, position, average luminosity 
and local metallicity will provide a definitive test of our pro- 
posal. 
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